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(N| Abstract 

The stochastic resonance (SR) in bistable systems has been extensively discussed with the use 






of phenomenological Langevin models. By using the microscopic, generalized Caldeira-Leggett 
(CL) model, we study in this paper, SR of an open bistable system coupled to a bath with a 



f^ nonlinear system-bath interaction. The adopted CL model yields the non-Markovian Langevin 

equation with nonlinear dissipation and state-dependent diffusion which preserve the fluctuation- 
dissipation relation (FDR). From numerical calculations, we find the following: (1) the spectral 

a 

power amplification (SPA) exhibits SR not only for a and b but also for r while the stationary 

probability distribution function is independent of them where a (b) denotes the magnitude of 

multiplicative (additive) noise and r expresses the relaxation time of colored noise; (2) the SPA for 

coexisting additive and multiplicative noises has a single-peak but two-peak structure as functions 

i i of a, b and/or r. These results (1) and (2) are qualitatively different from previous ones obtained 

VO by phenomenological Langevin models where the FDR is indefinite or not held. 
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I. INTRODUCTION 

In the last three decades since a seminal work of Benzi et. al. [I] , extensive studies have 
been made on the stochastic resonance (SR) (for a recent review on SR, see Refs. [21 E])- 
Theoretical studies on SR were made initially with the use of the overdamped Langevin 
model subjected to white noise. It has been shown that the signal-to-noise ratio (SNR) of 
nonlinear systems to a small periodic signal is resonantly enhanced under some optimum 
level of noise, whose effect is called SR [21 [3]. The resonant activation occurs when the 
frequency of an applied periodic signal is near the Kramer escape rate of the transition 
from one potential minimum to another. Subsequent studies on SR examined effects of 
colored noise (for a review on colored noise, see Ref. [3]). Calculations with the use of the 
universal colored-noise approximation (UCNA) [5] have shown that noise color suppresses 
SR monotonically with increasing its relaxation time r [SHE]. Ref. [9] has claimed that 
when t is increased, the SR is first suppressed for small r but enhanced for large r with the 
minimum at intermediate r. In recent years, studies on SR have been made by using more 
sophisticated noise: correlated additive white noise and multiplicative colored noise [TU1 - 
[13] • SNR of the Langevin model subjected to such noises is shown to exhibit a complicated 
behavior as functions of magnitudes of additive and multiplicative noises and the correlation 
strength between the two noises [rH - TH] . For example, SNR may have the bimodal structure 
when these parameters are varied. 

It should be noted that Langevin equations having been adopted for studies on SR are phe- 
nomenological ones which have some deficits such that a lack of the fluctuation-dissipation 
relation (FDR), related discussion being given in Sec. IV. We expect that in almost all real, 
non-isolated systems, the environment acts as a thermal bath which is a source of noise. 
The microscopic origin of additive noise has been proposed within a framework of Caldeira- 
Leggett (CL) system-bath Hamiltonian [ToT[TT] . The CL model with a linear system-bath 
coupling leads to the non-Markovian Langevin equation where memory kernel and noise 
term are expressed in terms of bath variables and they satisfy the FDR [T5l - fTT] . The FDR 
implies that dissipation and diffusion come from the same origin. By using the generalized 
CL model including a nonlinear system-bath coupling, we may obtain the non-Markovian 
Langevin equation with multiplicative noise which preserves the FDR [181 HH] • The nature 
of nonlinear dissipation and multiplicative noise has been recently explored with renewed 



interest [201422] . Nonlinear dissipation and multiplicative noise have been recognized as im- 
portant ingredients in several fields such as mesoscopic scale systems [2"31 |2"1] and ratchet 
problems J2SH2B]. 

Quite recently, we have studied the generalized CL model for a harmonic oscillator which 
includes nonlinear nonlocal dissipation and multiplicative diffusion [29|. It has been shown 
that the stationary probability distribution function (PDF) is independent of magnitudes 
of additive and multiplicative noises and of the relaxation time of colored noise [20"H2"2| 129] , 
although the response to applied input depends on noise parameters [29J . This is in contrast 
with previous studies [3"UH3"2"] which show that the stationary PDF of the Langevin model for 
the harmonic potential is Gaussian or non-Gaussian, depending on magnitudes of additive 
and multiplicative noises. We expect that when the generalized CL model with a nonlinear 
coupling [20H221 [29] is applied to a bistable system, properties of the calculated SR may be 
rather different from those having been derived by the phenomenological Langevin models 
[U EHH]. It is the purpose of the present paper to make a detailed study on SR of an open 
bistable system described by the microscopic, generalized CL model with nonlinear nonlocal 
dissipation and state-dependent diffusion |29j. Although the UCNA has been employed in 
many studies on SR for colored noise [TUHH] because it is exact both for r = 0.0 and r = 00 
[5], its validity within 0{r) is not justified [331 G2]- In this study, we investigate SR for a 
wide range of r value by simulations. 

The paper is organized as follows. In Sec. II, we briefly summarize the non-Markovian 
Langevin equation derived from the generalized CL model with OU colored noise [2"UH2"2"| |2"$] . 
In Sec. Ill, after studying the stationary PDF, we investigate SR of an open bistable 
system for three cases: additive noise, multiplicative noise, and coexisting additive and 
multiplicative noises. In Sec. IV, we examine a validity of the Markovian Langevin equations 
with the local approximation, which have been widely employed in previous studies on SR 
[U EHH]. Sec. V is devoted to our conclusion. 



II. THE GENERALIZED CL MODEL 



A. Non-Markovian Langevin equation 



We consider a system of a Brownian particle coupled to a bath consisting of iV-body 
uncoupled oscillators, which is described by the generalized CL model [2014221 129] , 



H = Hs + H b + Hi 



with 
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Here Hs, Hb and Hi express Hamiltonians of the system, bath and interaction, respectively; 
x, p and V(x) denote position, momentum and potential, respectively, of the system; q n , p n , 
m n and u n stand for position, momentum, mass and frequency, respectively, of bath; the 
system couples to the bath nonlinearly through a function 4>(x); f(t) expresses an applied 
external force. The original CL model adopts a linear system-bath coupling with 4>(x) = x 
in Eq. (J3J) which yields additive noise [17] . By using the standard procedure, we obtain the 
generalized Langevin equation given by [T7H20] [29] 



with 
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where j(t — t') denotes the non-local memory kernel and ((t) stands for noise, and dot and 
prime stand for derivatives with respect to time and argument, respectively. Dissipation 
and diffusion terms given by Eqs. (J5| and ([6]), respectively, satisfy the FDR, 



(C(t)C(O) = ksT^t-t') 



(7) 



where the bracket (-) stands for the average over initial states of g n (0) and p n (0) [18H20] . 



We adopt the OU process for the kernel j(t — t') given by 

l(t - f) = (^) e-l^l/-, (8) 

where 70 and r stand for the strength and relaxation time, respectively, of colored noise. 
The OU colored noise may be generated by the differential equation, 
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where £(£) expresses white noise with 

m) = o 5 (mw)) = 8(t-i!). (10) 

Equations rt9j) and (10) lead to the PDF and correlation of colored noise given by 

(11) 
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where /3 = l/k B T. 

In order to transform the non-Markovian Langevin equation given by Eq. (El) into more 
tractable multiple differential equations, we introduce a new variable u(t) [19], 122] . 

u{t) = - f 7 (t - t')<f)'{x{t'))x{t') dt', (13) 

Jo 

which yields four first-order differential equations for x(t), p(t), u(t) and £(£) given by 

x(t) = p(t), (14) 

p{t) = -V'(x) + <f>'(x(t)) u(t) + f(t) + <f>'(x(t)) C(t), (15) 

u(t) = -^-(^)<P\x{t))p{t\ (16) 
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From Eqs. (Jl4|)-(Jl7|), we obtain the multi-variate Fokker-Planck equation (FPE) for distri- 
bution of P(x, p, u, C, t) , 
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The stationary PDF with f(t) = is given by [19] 

P(x, p) = J P{x, p, u, C) du dC oc e -P[p 2 H+v(x)\ _ ( 19 ) 

B. Markovian Langevin equation 

It is worthwhile to examine the local limit of Eq. ([7]), 

(C(*)C(0)o = k BT lit - t') = 2k B T 70 S(t - 0, (20) 

which leads to the Markovian Langevin equation given by 

x(t) = -V / (x(t))- 7 o0 / (x(t)) 2 x(t) + v / 2^T0 / (x)e(t) + /(t), (21) 

preserving the FDR. It is evident that the Markovian Langevin equation becomes a good 
approximation of the non-Mar kovian one in the limit of r — > 0. 



From Eq. (21), we obtain first-order differential equations, 



x(t) = p(t), (22) 



p{t) = -V'(x) - l0 <f>'(x(t)) 2 p(t) + y/2k B Ti <j>\x{t)) £(t) + f{t). (23) 

The relevant FPE for the PDF of P(x,p,t) is expressed by 

dP ^ P,t) = -^ c pP(x,p,t) + ^ p [V'(x)-f(t)+i 0( f>'(x) 2 plP(x,p,t) 

+ k B T l0 — 4f(x)— <f>'(x)P(x,p,t). (24) 



The stationary distribution of Eq. (24) with f(t) = is given by 

P{x,p) oc e -^ 2 ' 2+v ^\ (25) 

which is consistent with the result of the non-Markovian Langevin equation given by Eq. 



(19). 



C. Overdamped Langevin equation 

We will derive the overdamped Langevin equation from the Markovian Langevin model 



given by Eq. (21). The overdamped limit of the Markovian Langevin equation is conven- 



tionally derived with setting x = in Eq. (21). This is, however, not the case because 



dissipation and diffusion terms are state dependent in Eq. (21). Sancho, San Miguel and 



Diirr [37] have developed an adiabatic elimination procedure to obtain an exact overdamped 



Langevin equation and its FPE. In order to adopt their method [37J, we rewrite Eq. (21 ) as 



x(t) = -V'(x(t))-\(x(t))x(t)+g(x(t))t(t) + f(t)i 
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By the adiabatic elimination in Eq. (26) after Ref. [37], the FPE in the Stratonovich 
interpretation is given by 
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where we employ the relation: g(x) 2 = 2ksT\(x) derived from Eqs. (27) and (28). The 
relevant Langevin equation is given by [37] 
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Note that the second term of Eq. (30) does not appear when we obtain the overdamped 



Langevin equation by simply setting x = in Eq. (26). It is easy to see that the stationary 



distribution of Eq. (29) with fit) = is given by 
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III. BISTABLE SYSTEMS 



A. Stationary PDF 

We apply the generalized CL model mentioned in the preceding section to a bistable 
system, where V{x) and 0(x) in Eqs. 0-((3]) are given by 



V{x) 

0(x) 
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ax 
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(32) 
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a and b denoting magnitudes of multiplicative and additive noises, respectively. The symmet- 



ric bistable potential given by Eq. (32 ) has two stable minima at x — ±1.0 and unstable max- 



imum at x = 0.0. The height of the potential barrier is given by A = V(0) — V(±l) = 1/4. 



First we show calculations of the stationary PDF. We have solved Eqs. (14)-(17) by 



using the Heun method [351 ES] with a time step of 0.001 for parameters of 70 = 1.0 and 



ksT = A = 0.25. Simulations of Eqs. (14)-(17) have been performed for < t < 1000.0 



with discarding initial results at t < 200.0. Calculated results to be reported are averaged 
over 10 000 runs with initial states of Gaussian-distributed x(0), p(0) and ((0) with (x(0)) = 
(p(0)) = (C(0)> = 0, (x(0) 2 ) = (p(0) 2 ) = k B T, and (C(0) 2 ) = k B Tj /r. An initial condition 
of w(0) = has been used in all simulations. 

Calculated marginal PDFs of P{x), P{p), P{u) and -P(C) for f{t) = are shown in Figs. 
[IJa), (b), (c) and (d), respectively, where four sets of parameters are adopted: (a,b,r) = 
(1, 0, 10.0) (solid curves), (1, 0, 0.1) (dashed curves), (0, 1, 10.0) (chain curves) and (0, 1, 0.1) 
(dotted curves). We note that P(x) and P(p) in Figs. [l](a) and (b) are independent of 
the parameters of (a, b, r) although P(u) and -P(C) in Figs. [IJc) and (d) depend on them. 
Calculated P(x) and P(p) are in good agreement with marginal PDFs obtained by 



P(x) = P(x,p)dpoce-P v{x \ (34) 

P{p) = f P(x,p)dxoce-P p2/2 , (35) 

which show that P(x) has the typical bimodal structure while P(p) is Gaussian. Equation 



(11 ) shows that P(() is the Gaussian PDF whose variance depends on r for fixed 70 and T. 
In contrast, Fig. [TVc) shows that P(u) is Gaussian PDF for additive noise but non-Gaussian 
PDF for multiplicative noise. Properties of stationary PDFs for bistable potential shown in 
Fig. 1 are similar to those for harmonic oscillators except for P(x) (see Fig. 1 of Ref. [29]). 



B. Stochastic resonance 



1. Spectral power amplification 



We investigate SR, applying a sinusoidal input given by 



f(t) = g cos ( — j = g cos w t, (36) 
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FIG. 1: (Color online) Marginal PDFs of (a) P(x), (b) P(p), (c) P(u) and (d) P(() for 
(a, 6, r) = (1,0,10.0) (solid curves), (1,0,0.1) (dashed curves), (0,1,10.0) (chain curves), and 
(0,1,0.1) (dotted curves) calculated by simulations with 70 = 1.0 and ksT = A = 0.25. P{x) 
and P(p) in (a) and (b) are indistinguishable among the four sets of (o, b, r). P(u) and P(C) for 
(a,b,r) = (1,0,0.1) and (0,1,0.1) are multiplied by a factor of five in (c) and (d) where PDFs are 
arbitrarily shifted for a clarity of figures. 



where g, T and uq (= 27r/T ) denote its magnitude, period and frequency, respectively. A 



sinusoidal input given by Eq. ( 36 ) yields an averaged output given by 



fi(t) = (x(t)), 



(37) 



where (•) denotes the average over trials. We evaluate the spectral power amplification 
(SPA) S defined by 
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FIG. 2: (Color online) (a) SPA (S) as a function of a period of an input (To) for (a, b, r) = 
(0.0, 0.1, 0.1) (solid curve) and (0.1, 0.0, 0.1) (dashed curve), (b) S as a function of To for (a, 6, r) = 
(0.0,0.5,5.0) (solid curve) and (0.5,0.0,5.0) (dashed curve). 

where f[u] and //[a;] are Fourier transformations of f(t) and //(£), respectively. Simulation 
procedures are the same as those for a calculation of stationary PDFs except for initial 
conditions of x(0) and p(0) which are taken to be x(0) = 1.0 and p(0) = 0.0. Direct 
simulations for Eqs. (Jl4|)-(Jl7|) have been made for sets of a, b, r, g and T averaged over 
10 000 runs. 

Figurel^a) shows SPA (S) as a function of a period of an external input (T ) with g = 0.05 
for two sets of parameters: (a,b,r) = (0.0,0.1,0.1) (solid curve) and (0.1,0.0,0.1) (dashed 
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curve), the former (latter) corresponding to additive (multiplicative) noise only. SPAs for 
both sets of parameters have maxima at T ~ 9. It is interesting that S for additive noise 
nearly coincides with that for multiplicative noise. This coincidence, however, is not realized 
for a larger r and/or b as shown in Fig. \2Xb) where similar plots of the To dependence of S 
are presented for different two sets of parameters: (a, b, r) = (0.0, 0.5, 5) (solid curve) and 
(0.5, 0.0, 5) (dashed curve). Resonances of S in Fig. |2^b) are again realized at about T ~ 9 
which coincides with the result in Fig. [2|a). This arises from the fact that an increase in the 
inverse of the Kramers rate with an increase of r from r = 0.1 to r = 5.0 is approximately 
compensated by its decrease with an increase of b from b = 0.1 to b = 0.5. Figures [21(a) and 
p?[b) clearly show the resonant behavior of S when a period of an external sinusoidal input 
To is varied. In the following, we will separately present simulation results for three cases: 
additive noise, multiplicative noise, and coexisting additive and multiplicative noises with 
fixed values of T — 10.0 and g = 0.05. 

2. Case of additive colored noise 

First we apply only additive colored noise (a = 0.0) to the system. Figure [31(a) shows 
SPA as a function of b for r = 0.1 (dashed curve), r = 5.0 (solid curve) and r = 10.0 
(dashed curve). With increasing r, the b value where S has the maximum is increased, and 
the maximum value of S (S max ) is gradually decreased with the wider distribution of S. 
Figure [3]^b) shows the r dependence of S calculated for b = 0.2 (dashed curve), b = 0.4 
(solid curve) and b = 0.6 (dashed curve). S max in Figs. [3|a) and[3|b) is observed at a larger 
b for a larger r and vice versa. 

Figure |4Fa) shows a schematic plot of S as functions of b and r, which is estimated from 
simulation results shown in Figs. |3ta) and [31(b). We note that with increasing b or r, a 
location of S max departs from the origin of (6, r) = (0.0,0.0) and the magnitude of S max is 
gradually decreased. Circles in Fig. [4|b) denote locations of S max in the b-r plane obtained 
by simulations shown in Figs. [3ta) and (b). S max nearly locates along the chain curve 
expressed by b = 0.075r + 0.1. Fig. W(c) will be explained shortly. 
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FIG. 3: (Color online) (a) SPA S for additive colored noise as a function of b with r = 0.1 (dashed 
curve), r = 5.0 (solid curve) and r = 10.0 (chain curve), (b) S as a function of r with b = 0.2 
(dashed curve), b = 0.4 (solid curve) and b = 0.6 (chain curve) (a = 0.0, g = 0.05 and To = 10.0). 

3. Case of multiplicative colored noise 



Figure [51(a) shows the a dependence of S when we apply only multiplicative colored noise 
with t = 0.1 (dashed curve), r = 5.0 (solid curve) and r = 10.0 (dashed curve). With 
increasing r, the position of S max moves to a larger a and S max is gradually decreased. 
Figure |5](b) shows the r dependence of S calculated for a = 0.2 (dashed curve), a = 0.6 
(solid curve) and b = 1.0 (dashed curve). 
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FIG. 4: (Color online) (a) A schematic plot of S as functions of b and r (a and r) for additive 
(multiplicative) noise, (b) Locations of S max for additive colored noise in the b-r plane obtained 
by simulations (circles), the chain line being expressed by b = 0.075r + 0.1. (c) Locations of S max 
for multiplicative colored noise in the a-T plane obtained by simulations (circles), the chain line 
being expressed by a = 0.09r + 0.2. 

We expect from simulation results in Figs. [5ta) and[5tb) that the a- and r-dependent PSA 
is given by Fig. ffla) where we read the b axis as the a axis. As in the case where additive 
colored noise is added, a location of S max departs from the origin of (a, r) = (0.0, 0.0) and 
Smax is gradually decreased with increasing a or r. Circles in Fig. [4jc) denote locations of 
S max in the a-T plane obtained by simulations shown in Figs. |5^a) and (b). S max nearly 
locates along the chain curve expressed by a = 0.09r + 0.2. 
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FIG. 5: (Color online) (a) SPA S for multiplicative colored noise as a function of a with r = 0.1 
(dashed curve), r = 5.0 (solid curve) and r = 10.0 (chain curve), (b) 5 as a function of r with 
a = 0.2 (dashed curve), a = 0.6 (solid curve) and a = 1.0 (chain curve) (b = 0.0, g = 0.05 and 
T = 10.0). 

4- Case of coexisting additive and multiplicative colored noises 



Next we study the case where both additive and multiplicative colored noises are applied 
to the system. Figure |6](a) shows the a dependence of S for r = 0.1 with b = 0.0 (dashed 
curves), b = 0.05 (dotted curve), b = 0.2 (chain curves) and b = 0.4 (solid curves). The 
dashed curve for b = 0.0 shows that S has the maximum at a ~ 0.14, as having been 
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FIG. 6: (Color online) SPA S for coexisting additive and multiplicative colored noises with (a) 
r = 0.1 and (b) r = 10.0: b = 0.0 (dashed curves), b = 0.05 (dotted curves), b = 0.2 (chain curves), 
b = 0.4 (solid curves), and b = 0.6 (bold solid curve) (g = 0.05 and To = 10.0). 

shown in Fig. pi a). With adding small additive noise of b = 0.05, S max is increased and its 
location moves to a smaller value of a ~ 0.05. With further increase of b (> 0.2), however, 
the magnitude of S is gradually decreased and its maximum becomes obscure. A similar 
a dependence of S for r = 10.0 is plotted in Fig. [6^b) where S is monotonously increased 
with increasing b. 

Figure [7^a) expresses a schematic plot of S as functions of a and b for the case of coexisting 
additive and multiplicative noises, which are estimated from results of simulations shown in 
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FIG. 7: (Color online) (a) A schematic plot of S as functions of a and b for coexisting additive 
and multiplicative noises, (b) Locations of S max in the a-b plane: circles and squares denote S max 
for r = 0.1 and r = 10.0, respectively, obtained by simulations, chain and double-chain curves 
being shown for a guide of the eye. Along dashed lines, a value in Fig. |6ja) or WLb) is changed for 
a given b. 



Figs. |6ja) and|6jV). S has a small value at the origin of (a, b) = (0.0,0.0). With departing 
from the origin by increases of a and/or b, S is first increased and then decreased after 
passing the maximum. Circles and squares in Fig. |7[b) denote locations of S max in the a-b 
plane obtained for r = 0.1 and r = 10.0, respectively, from Figs. |6](a) and[6^b). S max nearly 
locates along a chain (double-chain) curve for r = 0.1 (r = 10.0): in the case of r = 10.0 
locations of S max are scattered because maximum of S become obscure [Fig. [6^b)]. 

Dashed lines in Fig. fT^b) express traces along which a value is changed for a given b in 
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Figs. J6](a) and[6^b). With increasing b in Fig. [6^a) for r = 0.1, the magnitude of a-dependent 
S is decreased because a dashed line departs from the maximum expressed by a chain curve 
in Fig. [Wb), except for b = 0.05 for which the a-dependent S is increased. In contrast, with 
increasing b in Fig. mb) for r = 10.0, the magnitude of the a-dependent S is increased, 
particularly for a < 0.5, because a dashed line approaches the maximum expressed by a 
double-chain curve in Fig. [71(b). These explain the difference of the b dependences of S in 
Figs. 7(a) and 7(b). 

IV. DISCUSSION 

A. Underdamped Markovian approximation for additive noise 

We will examine the local approximation to the non-Markovian Langevin equation given 
by Eq. ^ for additive colored noise with (j)'(x) = 6, 

r t 

<-2 



x = -V'(x) - b 2 / j(t - t') x(t') dt' + b C(t) + /(£), (39) 

Jo 

where the nonlocal kernel and OU colored noise are given by 

(C(0C(O) = k B T>y{t - if) = (^) e- |4 - 4 ' l/T . (40) 



As mentioned in Sec. II B, when employing the local limit given by Eq. (20), 

(C(t)C(O) = k BT l{t - = 2k B T 7o 5(t - t'), (41) 

we obtain the Markovian Langevin equation given by 



x 



■V'{x) - 6 2 7o x + v^KW + /(<), (42) 



where £(t) denotes white noise given by Eq. (10). If an alternative local approximation 
given by 

7 (t - f) = 2 7o s(t - 1% mat)) = C^IjA e -\t-tf\/r f m 



is adopted, Eq. (39) reduces to the Markovian Langevin equation given by 

x = -V'(x) - 6 2 7o x + b C(t) + f(t). (44) 
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FIG. 8: (Color online) The r dependence of S for additive noise (a = 0.0) calculated by non- 



Markovian Langevin equation [Eq. (39)] (solid curves) and underdamped (UD) Markovian approx- 



imation [Eq. (44)] (dashed curves) with b = 0.2 (circles), b = 0.4 (squares) and b = 0.6 (triangles) 
{g = 0.05 and T = 10.0). 



We should note that the memory kernel and colored noise given by Eq. (43) do not satisfy 



the FDR for r ^ 0.0, while those given by Eqs. (40) and (41 ) hold the FDR. The Markovian 



Langevin equation given by Eq. (44) with OU colored noise [Eq. (43)] has been employed 
for a study on SR in Ref. [8]. 

Solid curves in Fig. 8 show the r dependence of S calculated by the non-Markovian 



Langevin equation [Eq. (39 )] for b = 0.2 (circles), 0.4 (squares) and 0.6 (triangle) with ksT = 
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0.25, 70 — 1-0, g = 0.05 and To = 10.0, which show the maxima as having been explained 
in Fig. 3(b). Dashed curves express relevant plots of S for (underdamped) Markovian 



Langevin equation given by Eq. (44) with OU noise given [Eq. (43)]. At r — >■ 0.0, results of 
the Markovian approximation agree with those of the non- Markovian Langevin equation as 
expected. With increasing r, S in the Markovian approximation monotonously decreased 
as shown in Ref. [8], although for b = 0.2, S is slightly increased at 0.1 < r < 0.4 but 
decreased at r > 0.4. It is evident that results of the Markovian Langevin equation given 



by Eq. (44) with the local approximation given by Eq. (43) are quite different from those 



of the non-Markovian Langevin equation given by Eq. (39) except for a very small r. 



B. Overdamped Markovian approximation for additive noise 



In order to gain an insight into SR of the overdamped limit of the non-Markovian Langevin 



equation given by Eq. (39) subjected to additive noise, we have performed simulations, 



employing a large value of 70 in Eq. (40). Figure 9 shows the r dependence of S for 70 = 1.0 



(solid curve), 10.0 (dashed curve) and 100.0 (chain curve) with a = 0.0, b = 0.2, g = 0.05 and 
T = 10.0. With a larger 7 , S(r) for r — > becomes smaller and the peak position for S max 
becomes larger. S for all the cases investigated has the maximum. Our calculation suggests 
that SR might be realized against r even in the overdamped limit of the non-Markovian 
Langevin equation. 

In previous studies (6J [7] , the overdamped Langevin equation given by 



x = -V'(x) + ((t) + f(t), 



(45) 



with the OU noise correlation. 



<C(*)C(0) 



D 



-\t-?\/T 



(46) 



has been extensively studied within the UCNA. Equation (45) may be obtainable from Eq. 



(44) with setting x = and an appropriate change of variables. 



Figure 10 shows the r dependence of S calculated by simulations for the overdamped 



Markovian Langevin equation given by Eq. (45) with OU colored noise [Eq. (46)] for 
D = 0.10 (chain curve), 0.25 (dashed curve) and 0.50 (solid curve). With increasing r, all 
calculated results show monotonous decrease, in consistent with the result of the UCNA in 
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FIG. 9: (Color online) The r dependence of S for additive noise calculated by the non-Mar kovian 



Langevin equation [Eq. (39)] for 70 = 1.0 (solid curve), 70 = 10.0 (dashed curve) and 70 = 100.0 
(chain curve) with a = 0.0, b = 0.2, g = 0.05 and Tq = 10.0. 



Refs. [HI U\- However, they exhibit no signs of the SR against r, which is in contrast with 
the result in Fig. [9] obtained by the non-Markovian Langevin equation for a large 70. 



The discrepancy between results of non-Markovian Langevin equation [Eq. (39)] and 



Markovian Langevin equations [Eqs. (44) and (45)] whose results are presented in Figs. 8, 



9 and 10, arises from the fact that the FDR is preserved in the former but in the latter. We 



should note that a variation of r in Eq. (40) means simultaneous changes of the relaxation 



time of colored noise and of the correlation time of a memory kernel because they are 
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FIG. 10: (Color online) The r dependence of S for additive noise calculated by the overdamped 



(OD) Markovian approximation [Eq. (45)] with D = 0.10 (chain curve), D = 0.25 (dashed curve) 
and D = 0.50 (solid curve) (see text). 



mutually related by the FDR. In contrast, a conventional local approach given by Eq. (43) 



or (46) includes a variation of r only in the colored-noise relaxation time with a vanishing 



memory. 



Ref. [9] has investigated the overdamped limit of Eq. (39), taking account of the memory 
in a kernel given by 



ksT^t - = (£(*)£(*')) = 2 7 o S(t - + (^) 



■\t-t'\/r 



(47) 



where 7 and 71 express magnitudes of local and non-local dissipations, respectively. Equa- 
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tions (39) and (47) yield the over-damped Markovian Langevin equation given by 

,-t 



lox 



.y /(x) _ Tix + 7i*(0)e-«/- + fn 



^- t,)lT x{t')dt' + i{t) + f{t). (48) 



It has been shown that when r is increased, the SR is first suppressed for small r but 
enhanced for large r with the minimum at intermediate r [9J. Although the FDR is held in 



Eq. (47), its expression cannot be justified from a microscopic point of view. 



C. Overdamped Markovian approximation for additive and multiplicative noises 



For additive and multiplicative noises with 4>'(x) = ax + b, Eqs. (27), (28) and (30) in 
Sec. II C lead to the overdamped Markovian Langevin equation given by 

V'(x) k B Ta I 2k B T 



x 



+ 



7o(ax + 6) 2 7o(ax + 6) 3 y 70 (ax + b) 



t(t), 



(49) 



where £(t) is white noise given by Eq. (|10|). Equation (49) is, however, quite different from 



a widely-adopted phenomenological Langevin model given by 

x = -v'(x)+£(t)+xfj(t), 



(50) 



with correlations, 



(£(*)£(*')) = 2D 5(t - if), (f)(t)v(t')) = 2An 6(t - t') 

mm) = o, 



(51) 



where D and D m stand for magnitudes of additive and multiplicative noises, respectively. 



The stationary PDF obtained from the PFE for Eqs. (50) and (51) in the Stratonovich sense 
is given by 



InP(x) 



_jm_ dx -(± ]HD+DmA . 



(52) 



For the parabolic potential of V(x) = x 2 /2, Eq. (52) yields Gaussian or non-Gaussian PDF, 
depending on D and D m [3"0TI3"2"] . For the bistable potential of V(x) = x 4 /4 — x 2 /2, Eq. (52) 
leads to 



P(x) oc [D + D m x 

-(x 4 /<±-x 2 /2)/D 



oc e 



2,-UH t -D m -D)/2I%, e -a*/2D m for jy _, Q ^ ^ ^ Q ^ 

for D m = 0.0, 



oc \ x \VDm-i e -x 2 /™ m for ,0 = 0.0. 



(53) 
(54) 
(55) 
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Thus the stationary PDF of Eq. (52) depends on magnitudes of additive and multiplicative 



noises, in contrast with that of Eq. (31) which is independent of them. 

Studies on SR have been made in Refs. [TUrdl], by using the Langevin equation given 



by Eq. (50) but with generalized correlations, 



mm) 
mm) 



2D8(t-t'), (f}(t)f}(t')) = (D m /r) e 



-\t-t'\/r 



2Xy/DD m 6(t-lf), 



(56) 



where A is the coupling strength between additive and multiplicative noises. SNR calculated 
with the use of the UCNA and the functional integral method shows a complicated behavior 
as functions of D, D m and A. However, the phenomenological Langevin equation given by 



Eq. (50) with Eq. (56) has no microscopic bases and does not agree with microscopically 



obtained expression given by Eq. (49). 



V. CONCLUDING REMARKS 



Phenomenological Langevin equations mentioned in the preceding section have following 
deficits: 
(a) the FDR does not hold between a local dissipation and nonlocal colored noise in Eq. 



(43) for the Langevin equation given by Eqs. (44) or (45) except for r = 



(b) the nonlocal kernel and colored noise in Eq. (47) for the Langevin equation given by 



Eq. (48) have no microscopic bases although they hold the FDR, and 



(c) the FDR in the Langevin equation given by Eq. (50) with additive and multiple noises 



whose correlations are given by Eq. (51) or (56) is not definite 



The non-Markovian Langevin equation given by Eq. Q derived in this study from the 
generalized CL model is free from the deficits (a)-(c) in phenomenological Langevin models. 
We have studied the open bistable system by simulations, obtaining the following results: 
(i) calculated marginal PDFs for x and p are given by P(x) oc e~^ v( - x ^ and P(p) oc e~^ p / 2 , 
respectively, independently of noise parameters of a, b and r (Fig. 1) although SPA depends 
on them, 

(ii) SPA exhibits SR for a variation of r [Figs. |3tb) and[5tb)] besides those of a and b [Figs. 
|3]^a) and[5](a)], 
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(iii) SPA for coexisting additive and multiplicative noises shows unimodal but bimodal 
structures as functions of a, b and/or r [Figs. HJa) andfTta)], and 



(iv) the local approximation [Eq.(43)] to non-Markovian Langevin equation cannot provide 
a reliable description except for a vanishingly small r (Fig. pi). 

Items (i)-(iii) are in contrast with the results obtained by phenomenological Langevin models 
[SHU] where the FDR is indefinite or not preserved. In particular, the r dependence of SR 
in the item (ii) is different even qualitatively from previous results: a monotonous decrease 
of SNR with r JHHH] and the minimum of SNR at intermediate r [9]. It would be interesting 
to apply the present approach to a study of relevant phenomena in bistable systems such as 
the first passage time and resonant activation. 
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